YITP-SB-10-11 

Computation of the Soft Anomalous Dimension Matrix in Coordinate Space 

Alexander Mitov, George Sterman, Ilmo Sung 
C.N. Yang Institute for Theoretical Physics, Stony Brook University, Stony Brook, New York 11794-3840, USA 

(Dated: May 26, 2010) 

We complete the coordinate space calculation of the three-parton correlation in the two-loop 
massive soft anomalous dimension matrix. The full answer agrees with the result found previously 
by a different approach. The coordinate space treatment of renormalized two-loop gluon exchange 
diagrams exhibits their color symmetries in a transparent fashion. We compare coordinate space 
calculations of the soft anomalous dimension matrix with massive and massless eikonal lines and 
examine its nonuniform limit at absolute threshold. 

I. INTRODUCTION 

o 

As a result of recent investigations, our knowledge of the infrared (IR) singularities of massive gauge amplitudes has 
^^\ progressed from one [ll, 0] to two loops [j-Sl • We are now in a position to predict the single IR poles of any two-loop 
amplitude with arbitrary numbers of partons of arbitrary masses. This helps determine the behavior of amplitudes 
close to kinematic thresholds, and we now have available the two-loop input necessary for threshold resummation to 
next-to-next-to leading logarithm for heavy quark production [7|-|9[ . Much of the new information is encoded in the 
fvj ' two-loop anomalous dimension matrices for coupled massive and/or massless partons, which we refer to collectively 
as the "massive soft anomalous dimension matrix" [3|-[8|- 

A goal of the present paper is to clarify the relationship between two calculations of a key component in the two- 
loop massive anomalous dimension matrix, the color antisymmetric three-parton correlation. Implicit expressions 
7*^1 were presented for these correlations in [4, partly in terms of integrals in Euclidean space. Subsequently an elegant 
Oh analytic expression was presented in Ref. [a]. In Ref. I0| calculations were carried out in momentum space, and Ref. [J| 
^5 ■ in position space. Of course, these two approaches should give equivalent results, and we show below that when all 
contributions are taken in account, they indeed do. Part of our motivation in presenting the details leading to this 
expected outcome is that an apparent discrepancy in the calculations was raised in Q. This involves the diagrams 
called "double exchange" in Ref. [j| and "planar" in [6], illustrated here in Fig. [TJ As we shall see, this apparent 
^ \ disagreement arises simply because Ref. [4| exhibited results for these diagrams before renormalization. We present 
the remaining analysis here not as new results, but to confirm the equivalence of the two calculations, and because 
results of Ref. [4| were used in deriving the expression for the total cross section given in [^ . Of particular relevance 
^P . was the conclusion that the two-loop massive anomalous dimension matrix is diagonal in the s-channel singlet-octet 
basis for pair production from incoming light quarks or gluons at ninety degrees in the partonic center of mass. 

We will show in Sec. |ll] that apparent discrepancies between results presented in Refs. [J| and^y| are entirely due 
to two-loop diagrams with one-loop counterterms. These diagrams, which were not presented in 4], do not affect the 
color diagonalization of the anomalous dimension matrix at ninety degrees. This diagonalization was indeed confirmed 
in [6|. It is of interest to see the simplicity of the coordinate space analysis, which may have applications at higher 
orders [IMl- 

In Sec. mil we discuss similarities and differences encountered in the coordinate space calculations when all lines 
5—1 ■ are massless. Here, we present a discussion in which both one- and two-loop massless integrals are regularized 
dimensionally, and reproduce from a purely coordinate space analysis the absence of three-eikonal color correlations, 
a result found in [lj|, which relied on arguments in both momentum and coordinate space. 

Finally, in Sec. IIVI we discuss the limit of absolute threshold (pair creation at rest). Reference [6] discovered a non- 
uniform limit of the anomalous dimension matrix at absolute threshold, when expressed in terms of the scattering 
angle. Here we rederive the result found Ref. Q, using the methods introduced in Ref. [41 . We show how the 
singularity structure of the integrals restricts the non-uniform limit to a region of momentum space where the eikonal 
approximation fails. This limit of the soft anomalous dimension matrix is not directly relevant to the total pair 
production cross section [7|, |8|. For completeness, explicit expressions for the relevant integrals are exhibited in an 
appendix. 

II. TWO-LOOP DOUBLE EXCHANGE DIAGRAMS 

The formalism for the computation of soft anomalous dimension matrices has been discussed extensively in Refs. [l| 
and elsewhere. The relevant diagrams consist of radiative corrections to incoming and outgoing partonic propagators in 
eikonal approximation. The resulting diagrammatic expressions are scaleless, and vanish in dimensional regularization. 
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FIG. 1: Double exchange diagrams discussed in the text. The shaded circle represents an n-eikonal vertex. 




FIG. 2: Non-planar diagram for the two-loop anomalous dimension matrix. 

Their ultraviolet poles, however, define counterterms and hence the anomalous dimensions we are after |ll.ll3|. In this 
formalism, we factorize the full eikonal amplitude into a soft function, with one logarithmic singularity per loop and 
all nontrivial color exchange, and color-diagonal jet functions. This factorization ensures that in the limit of vanishing 
eikonal masses the soft function is free of collinear singularities [l3|, [ij] . A natural generalization of this procedure 
to eikonals of arbitrary mass is the scheme employed in Ref. |al, i n which the jet functions are defined as the square 
roots of the low-mass limits of the corresponding form factors |15l.[l6l|.^ 

The most challenging two-loop contribution to the massive anomalous dimension matrix is the set of non-planar 
diagrams with a three-gluon vertex connecting three eikonals, of the type shown in Fig. [5] Although these diagrams 
are rather complex, they have no UV-divergent subdiagrams. It was shown in Ref. [l3| by a momentum-space change 
of variables, that the ultraviolet poles of these diagrams vanish when all three of the interacting eikonals are massless. 
This result was rederived in Ref. [J| in coordinate space and extended to diagrams with two massless interacting 
eikonals, while a manifestly nonvanishing integral representation was given for the UV pole when all three interacting 
eikonals are massive. An analytic expression for the "non-planar" diagrams like Fig. [2] with massive external partons 
was derived in Ref. [6]. Surprisingly, it turned out to be very simple (see appendix). In numerical tests, it agrees with 
the integral representation given in Ref. [j| . 

In this section, we revisit the double-exchange diagrams of Fig.[T]in coordinate space, including their renormalization. 
As explained in Ref. [J| and confirmed below, the contributions from the two-loop double exchange diagrams are fully 
symmetric in color after combining pairs like those shown in Fig. [TJ They contain, however, subdiagrams with UV 
poles, and thus require one- loop counterterms. These counterterms were not considered explicitly in Ref. [j|. In Fig. [3] 
we show the renormalization of the double exchange diagrams of Fig. [l] with divergent loops replaced by counterterms. 
In Fig. [5^, for example, the counterterm depends on the pair of momenta Pk,Pj and the remaining one- loop diagram 
on Pi and Pj . 

To determine the counterterms in our coordinate space formalism, and to establish notation, we review the one- 
loop case, illustrated jxi Fig. |4l The counterterms, in turn, are found from the poles of one-loop diagrams. Following 
the notation of Ref. ^ (see also Appendix B of Ref. [Ijl for the same calculation in slightly different notation), we 



^ We note that Ref. '9| analyzes the massive amplitude direetly, without subtractions, but direet eomparison ean still be made, because 
by construction subtractions do not affect color exchange. 
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FIG. 3: Order g'* diagrams containing one-loop UV counterterm. 




FIG. 4: A gluon exchanged between an incoming quark i and an outgoing quark j that contributes to Eq. ([2]l. 



introduce the scalar propagator in coordinate space, 

A{x-y) = i /"d4-2e^g-»fc.(x-a) 
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where the third relation defines the constant c(e) = I + 0{e). 

The coefficient of ctg/T^ in the one-loop correction corresponding to Fig. |4]is a matrix in color space, ^ 
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The color factors associated with the gluon exchange are represented in the basis-independent notation introduced in 
Ref. [l3| and generally employed in analyses of soft anomalous dimensions. The final relation in Eq. ([2]) shows the 
contribution of the UV pole to the one-loop anomalous dimension matrix in terms of the function 
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^ This factor is referred to as the "velocity factor" in Ref. [li|. Because all partonic lines are eikonal, we need not distinguish between 
velocities and momenta for the external lines. 



with mjj = p^j ^ 0. The massless hmit is explored in Sec. IIIII For the purposes of this discussion, we may take 
—•pj ■ Pi > 0, corresponding to one hne incoming and the other outgoing. In this case, I{pi,pj,e) is reaL The exphcit 
results for the functions I^^^'^'{pi,pj) for both signs of the invariant are given in the appendix. 

It is important to note that the prefactor in Eq. ([2|) is arbitrary at 0{e), because it reflects the choice we make 
to define the integral over A-,-, a variable with dimensions of length squared. In effect, our choice for these functions 
constitutes a minimal subtraction scheme. In particular, the factor nij^ in ([2]) and following expressions is kept as a 
bookkeeping device only, because it matches the definition of I{pi,pj,e). For the purposes of calculating anomalous 
dimensions it is only necessary that the prefactors for all such integrals be defined in the same way. Then, as we will 
see below, overall factors like c(e) and m^"^ will not affect our results for the dimensionless pole terms of A^*^^^ 

The one-loop soft anomalous dimension matrix is defined only up to color-diagonal contributions associated with 
its eikonal lines, which absorb coUinear logarithms and are factorized from the soft function. In any such factorization 
scheme, the one-loop anomalous dimension matrix is found from the residues of the poles l/(— 2e) in the one-loop 
amplitudes of the soft function [1^. As noted above, we use the scheme of Ref. [^, in which the soft function is the 
full eikonal amplitude divided by the low-mass limit of the square root of the form factor for each external eikonal 
line. The one-loop soft anomalous dimension matrices found from the A^,-,, are then given by 
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where in the first line we subtract for each external eikonal a term that is collinear-singular in the massless limit, and 

where Cq = Cq = C, 

logarithms, as may be readily verified from the expressions for /' 

of Eq. dU, we recall the color identities, [13] 



re Uq = Uq = Up, Cg = Ca refer to the color content of the eikonal lines. This procedure eliminates all collinear 



given in the appendix. To verify the second line 
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with i = F^A. We are now ready for an analysis of the double-exchange diagrams. 

The coefficients of (a^/Tr)^ in the two-loop exchange diagrams of Fig. [T^, involving lines i, j and k, can be written 



in terms of the coordinate space integrals, /, Eq. ([3]), as 
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where the variables Aj and A' represent the locations of the vertices along eikonal line j, to which two gluons are 
attached, as in Fig. [T] The ordered pairs of arguments Pi,Pj {pj,pk) denote the eikonal lines connected by the outer 
(inner) gluon exchange of the diagram. The double pole of Eq. (jH]) is associated with the independent limits Xj —?■ 
and X'j — )■ of distances from the origin to the vertices along the j eikonal. The factors I{pi,pj,e) and I{pk,Pj,e) are 
the result of integrals like Eq. ([3]), which are over the ratios of origin-to- vertex distances. These integrals converge at 
both zero and infinity for fixed masses and momentum invariants. Thus, the individual poles arise from independent 
scaling of the two vertices of the ij- or /cj-exchange in Fig. [1] to the origin. Any remaining single pole represents a 
scaling of the four vertices of the entire diagram together to the origin. The latter defines the diagram's contribution 
to the two-loop anomalous dimension. 

The UV poles of the Xj integrals in Eq. (jS]) can readily be isolated, because the functions I{pi,pj,e), defined in 
Eq. ([3]), are independent of the A^'s. We find. 
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with both double and single poles, in a form that is clearly symmetric in the pairs ij and kj [J|. 

The complete set of two- loop diagrams corresponding to Fig. [1] includes as well counterterms for the UV-divergent 
subdiagram consisting of the (inner) exchange as shown in Fig. [31 For Fig. [3^, this one-loop counterterm is simply 



—A4^^\pk,Pj,s), defined by Eq. ([2]). For convenience, we include the factor c(e) as part of the renormalization scheme. 
The diagrammatic representation of the corresponding two-loop contribution is shown in Fig. [5^, and is given by 
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again providing both single- and double-pole contributions. Combining Fig. [T^ with its one-loop counterterm, Fig. 
[5^, we find 
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In this expression, we see directly that the double pole is symmetric in the momentum pairs pi , pj and pk , Pj , while 
the single pole is purely antisymmetric. 

We next combine the foregoing result with the diagram in which the roles of the pairs i] and jfc are exchanged, as 
in Figs.[T|D andEh). Including color factors, we find the following pole and color structure, 

M<^'\p.,Pj,Pk,e) + M^^^^''-\p.,Pj,Pk,e) + (z ^ k) = - {T, • T„T, • TJ ''^ I^-'\p,,p,) I^-'\pk,Pj) 
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We recognize that the single pole terms in this expression enter entirely with antisymmetric color factors, while the 
double poles are entirely symmetric in color [j|. 

The combination of all such three-eikonal double exchanges (DE) can now be written down for an arbitrary number 
of eikonal lines. The double pole terms of Eq. ([T0|) are cancelled by double poles of two-loop counterterms, which we 
see are the expansion of exponentiated one-loop counterterms. This leaves the single-pole, color-antisymmetric terms. 
Recalling |l3j that the contributions of these diagrams to the two-loop anomalous dimension matrix is the residue of 
the single pole, — 1/(4£), we find 



^(3£;,DE) _ 



(t) 5 E •/ 



i>j>k=l 



abcrj.arj.brj.c ^ ^ 

I,J,K={i,j,k) 



UK 



I(-'Hpi,Pj)I'^"Hpk,Pj), 



(11) 



where we have exhibited the antisymmetric color structure explicitly by carrying out the commutators in Eq. (jlOp . 
In this sum, the e-expansions of the overall factors Eq. ([TU)) . c^{e)m'j'^ , cancel because of the antisymmetry in (|lll) . 

which is invariant under any modification I^'^'{pi,pj) — > I^^'{pi,pj) + const x I^~^'{pi,pj). In fact, we can generalize 
this result to any overall rescaling of the integrals /, 
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where C has an expansion, C{pi,pj,e) = 1 + eC'-^\pi,pj) + As can be seen immediately from Eq. ([TU)). any such 

overall rescaling has no effect on the antisymmetric single pole terms of the two-loop anomalous dimension. Such a 
rescaling, of course, does change the symmetric color terms at the single-pole level, and should be thought of as a 
change in scheme for the soft function. Equation (ITTT) is in the form found in Ref. fo] . The exact expressions for the 
integrals I(pi,pj, e), in terms of which our results have been expressed so far, are given explicitly in the appendix for 
both space- and time-like kinematics. These results confirm the consistency of the two approaches in Refs. [J and Q 
that we set out to establish. 

To conclude this section, we note for completeness the form of the full three-eikonal anomalous dimension matrix, 
which is found by combining the results from the double-exchange diagrams, Eq. dill) , with the three-gluon diagram. 
Fig. [2j The latter has the same color structure, and can be incorporated into Eq. (fTTj) by simply adding terms to the 
factors /(°) iuEq. ^ 3> 
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The explicit result for I^g can be found in the appendix. Finally, we note that to construct the full two-loop anomalous 
dimension matrix for an n-point amplitude one needs to combine Eq. (J13p with the diagrams in which gluons connect 
to only two eikonals. The latter, which have the same color structure as the one-loop matrix, have been given explicitly 
inRef. [|. 

III. MASSLESS EIKONALS IN COORDINATE SPACE 

The diagrams treated in the previous section, of course, also appear in the calculation of the anomalous dimension 
matrix with massless eikonal lines [1^. In Ref. [l^, the same double exchange diagrams considered above, Figs. 
[T^,b, were shown for the massless case to give vanishing contributions to the two-loop anomalous dimension matrix 
using a momentum-space analysis. That is, their only contributions are cancelled by two-loop counterterms that 
are expansions of exponentiated one-loop counterterms. It is worthwhile examining how this result is maintained in 
coordinate space when taking into account the counterterm diagrams of Fig. |31 

We begin by recalling that the soft function for massless lines must be defined to eliminate coUinear poles on a 



loop-by- loop basis. We choose the form employed in Ref. 13|, and mentioned above, where the soft function is defined 



as the ratio of the full eikonal amplitude to the product of the square roots of color-singlet form factors, one for 
each eikonal line. This is the analog of the color diagonal subtraction in Eq. ^ in the massive case. We will not 
need the details of this construction here, but observe that it requires that at each loop order, double poles appear 
with only color-diagonal coefRcients after the sum over diagrams. We will see this feature emerge explicitly below, 
as we discuss the color structure of one- and two- loop exchange diagrams, beginning with the massless analog of the 
one- loop diagram, Eq. ([2]). 

For computations involving massless eikonals, we continue to use momenta pi in place of dimensionless four- 
velocities. Because the integrands that define the amplitudes are scaleless in the velocities, this is a trivial substitution 



before integration. For pf = p"^ = we then have 
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where we have multiplied by /i^^ so that both sides of this expression are dimensionless. This will help motivate our 
choice of renormalization scheme for this integral. As in our discussion of the massive case, we treat —2pi ■ pj as a 
positive quantity, and then analytically continue to negative values. The coordinate space expression in p^ has two 
scaleless integrals rather than one, as in Eq. ([2]). The treatment of these integrals requires some discussion, even in 
the one-loop case. 

The presence of mixed infrared and collinear singularities in a Lorentz-invariant integral make it impossible to 
maintain at the same time the correct power behavior under the scaling of momenta and invariance under rescalings 
of the velocities, a feature sometimes referred to as the "cusp anomaly" [12]. As a result, the prescription to identify 
the ultraviolet pole does not, as in the massive case, uniquely fix the coefficient of that pole. We can organize these 
ambiguities by changing variables from the A's, which have dimensions of length squared, to dimensionless variables. 
This is a unique change of variables in this case, if we insist that the integral (fTl)l depend only on pi ■ pj . We are 
therefore led to the dimensionless variables 

h = ^3 {-'^Pi ■ Pj) 1 k^h {-'2pi ■ Pj) ■ (15) 

After this change of variables, the dimensional content of the integral is manifest. To isolate the ultraviolet pole, we 
use the symmetry between li and Ij, and write 
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The li integral is carried out at fixed Ij for infrared regularization (e < 0), to give a simple pole. Next, the ultraviolet 
pole of the Ij integral is isolated. This results, of course, in an expression with both single and double poles. 
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Momentum-dependence appears only in single poles, while double poles are momentum- independent. In either form 
of ([T7|. we can revert to a velocity dependence by choosing j5i = Pi/n- Keeping the momentum-dependence, we can 
identify the first line of Eq. (fT7| as the eikonal contribution to the familiar function I*^^' defined in [13] ■ As in the 
massive case above, the choice between the first and second lines is a choice of renormalization scheme. In fact, we 
could use the first line to define an alternative scheme to define the soft function. In this discussion, however, we stick 
with the choice of pure poles, as in Ref. |13|. 

As in the example of Eq. ((31), color conservation, Eq. (IS|), implies a form in which, as anticipated, the double poles 
organize themselves into color-diagonal terms. 
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The order as double pole terms cancel against the form factors in the definition of the soft function. The resulting 
one-loop counterterms that correspond to Eq. p^ are then pure single poles. 

Going on to two loops, wc follow the reasoning leading to Eq. (ITTT) in the massive case. The essential difference is 
that all integrals remain scaleless, and we do not encounter convergent integrals like Eq. ([3]), whose expansion in e 
can give single poles at two loops. For the massless case, the two-loop integrals produce no expansion in e beyond 
the pole terms that define the scheme. Alternately, there is no pole associated with a single scaling of all the vertices 
in the diagram together, of the sort we identified in Eq. ^ with massive eikonals. That is, in the massless case each 
integral of the two-loop diagram is fully independent. The massless analog of Eq. ([6]) for Fig.[T|D is then. 
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where we assume all lightlike momenta, and where the ordering of the arguments of Ai '^^ determines the ordering of 
the color matrices. The integrand here is manifestly symmetric in pi and pk- As a result. Fig. [1^ differs from Fig.[T|D 
only by the orderings of their integration parameters Xj and A'- and of their color factors. A simple exchange of the 
two integration labels eliminates the antisymmetric part in the sum of the two diagrams, and we find 
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Only the symmetric part survives in the sum of the two diagrams, considered as color matrices. This is the coordinate 
space analog of a momentum space argument leading to the same result in Ref. [l3|- As in the massive case, we 
are left to analyze potentially antisymmetric contributions from diagrams like Figs. |3l with counterterms. Unlike the 
massive case, however, the one-loop integrals of these diagrams are scaleless in both integration variables, and there 
are no finite terms, beyond those that may have been included in the scheme that defines the one-loop counterterms. 
This ensures that the counterterm diagrams are proportional to the original diagrams, and that their sum remains 
symmetric. In principle, these one-loop counterterms include double as well as single poles, but as we have seen in our 
one-loop example, double poles systematically cancel in the soft anomalous dimension matrix, a result that extends 
to ah orders ^^. 

The origin of the difference between the massive and massless cases is easy to trace to the one-loop integrals of Eqs. 
^ and (|16p . respectively. The integrands in the two cases differ only when the scaled integration variable a = Xi/Xj 
in ([3]) is large (or small) enough so that 
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Only for these "collinear" regions do the integrands distinguish between massless and massive eikonal lines. For any 
fixed mass, however, these regions are always present. Within each diagram, therefore, the limit of zero mass does 
not commute with the integrations. This is not the case for the full soft function, however, because these collinear 
re gion s, whether regularized by dimensions or masses, cancel in the ratio of the full eikonal amplitude to jet functions 
[l|, [la, [201 , after the sum over a gauge-invariant set of diagrams. 





FIG. 5: Two-loop double-exchange diagrams and their one-loop counterterms at order (j* in the two-to-two scattering reactions 
qq — s- QQ and gg — >■ QQ. 



IV. ANOMALOUS DIMENSIONS NEAR ABSOLUTE THRESHOLD 



In this section, we discuss the interpretation of threshold behavior of the eikonal soft function representing a two- 
to-two reaction of the type qq — > QQ or gg — >■ QQ, as in Fig. \5\ We will concentrate on the effects of the double 
exchange diagrams shown there, involving (for example) an incoming light quark or gluon, labelled 1, and the two 
outgoing particles, labelled 3 and 4. Our goal here is to rederive and interpret the non-uniform behavior in the massive 
anomalous dimension matrix found in Ref. 6] at absolute threshold, where s — > 4M^, with Af the heavy quark mass. 
We will confirm and interpret the conclusion 8] that this behavior relates to a region in momentum space in which 
the eikonal approximation does not apply to the underlying partonic cross section. 

Equation (llip above shows how the double-exchange diagrams contribute to the anomalous dimension matrix. For 



applications to heavy quark production, as illustrated in Fig. [5l we take pi — P2 — 
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PI = M^. We consider 



the incoming mass m to be small, m << M, and with heavy quark production in mind will eventually take it to zero. 
To find the threshold, s — {pi +^2)^ -^ 4M^, behavior for this two-to-two process, we take the /3 — \J\ — ^Ivp / s — ?> 0, 
with P the center of mass velocity for the heavy pair. ^ Specifically, using the expressions given in the appendix for 
in -^ M we find for the terms in Eq. (fTTj) with i = 4:, j — 3, k — 1 corresponding to Fig. \5\ 



l'^'\p^,p,)l^-'\p,,p,)-l'^-'\p^,ps)l^"\pup,) 



In 



ITT 



n f niM\ n^ 
6 



\—) 



i7rln(4/3) 



0(/?), (22) 



where ti = — 2(pi ■ p^) and ui = — 2(pi -^4) are customary invariants. 

As noted in Sec. IIIII above, the double logarithm in Eq. (|22p . which originates from the collinear region, must be 
cancelled by the diagram with a three-gluon interaction ^], which we reproduce in the appendix, Eq. (jA12p . After 
this cancellation, in the threshold limit the diagrams of Fig. [5] then give a contribution to the anomalous dimension 
matrix. 



^(431,flg5) 



1 



, ^^afccrparpbrj, 



In 



fmM 



i7rln(4/3) 



ITT \ TT^ 



0{p) 



(23) 



^ In terms of the formulas given in the appendix, this corresponds to ^34 
(equivalently, 1113 — >■ 1) in Eqs. | |A3| | and | |A4| |. 



1 (equivalently, 1134 — )■ 0) in Eqs. l)A9|l and jAllll and xi 



This result is still infrared divergent in the limit m — >■ 0. When we add the diagrams in which the roles of lines 3 and 
4 are reversed, however, we find the same infrared-finite Tg given in Eqs. (53)-(55) of Ref. [6j], 



.(431) ^ 1 



r^*-"'^ = - ir'"=T1TlTl In ( -^ 



1 / TT^ 

- (^-- - z^ln(4/3) 



OW) . (24) 



Equation ((24)) has a very interesting behavior in the limit /3 ^ 0. By adopting the center of mass relation, ti = 
— s(l — /3cos9)/2 and similarly for ui, we observe that 



rr" = - ir'^TlT'.T'i cos9 [^ + z^ln(4/3) ) +0(/3) , (25) 



in which the Coulomb singularity combines with an angular-dependent factor that vanishes at threshold to give a 
term that remains finite for /3 ^ 0. As pointed out in Ref. [6|, this term depends on the center of mass scattering 
angle, no matter how small (3 is. The /? — ?► limit would thus appear to be ambiguous. What are we to make of this? 

The factorization formalism that leads to the soft anomalous dimension matrix [l|, |20| applies in the approximation 
that the momentum of the active partons is much larger than that of the soft radiation which the anomalous dimension 
matrix is used to resum. The use of an anomalous dimension matrix thus requires that the eikonal approximation 
apply, that is, that the emission or absorption of this radiation leaves the four- velocities of the heavy quarks essentially 
unchanged. In fact, in the momentum region that produces the Coulomb singularity, the center of mass kinetic energy 
of the quark pair is order to/3^ in virtual as well as real states. Soft radiation emitted by the produced pair can 
carry an energy of no more than this order without violating the eikonal approximation and washing out the 1//3 
dependence. Thus, as /? vanishes, the energy range of virtual or real radiation to which we can apply the eikonal 
approximation vanishes even faster, as the square of the relative velocity. Although the /3-independent piece of the 
anomalous dimension in (|25p . which couples singlet and octet color states, is present at any finite value of f3, its range 
of applicability shrinks to zero for /3 — )• 0. Correspondingly, as argued in Ref. Q, in the inclusive cross section, as 
opposed to an elastic amplitude, it is only the range of energies m > uj > mffl that contribute to threshold logarithms. 

The Coulomb singularity at zero relative velocity /3 is a characteristic feature of heavy particle production and the 
analysis of bound-state formation [21|. We can trace the origin of the 1//3 dependence in a particularly simple fashion 
for the integrals I{pi,pj,e), defined in Eq. Q for the exchange diagrams. Setting for dcfiniteness j = 3,fc = 4, the 
integration variable a = A3/A4 measures the relative distances from the origin along the two heavy-particle eikonal 
lines, with momenta ps and p^. The a integration contour encounters two singularities, at 



P3-P4 , 1 /7 ^2 4 , . 

= l±2p±ie +0{P^) . (26) 



In the limit of vanishing relative velocity, then, the a integral is pinched between coalescing singularities at cr = 1, 
whose separation vanishes linearly with j3. In terms of the original integrals over A3 and A4, the singularity at 
cr = 1 ± 2/3 corresponds to singularities in (say) A3 at A3 = A4(l ± 2/3). In particular, for /3 = 0, the two ends of 
the gluon propagator are at the same point in coordinate space, X3P3 = A4P4. We can use this analysis to interpret 
further the range of applicability of the anomalous dimension matrix. 

Working in the P3, p4 center of mass frame, the quantities A3 + A4 and A3 — A4 are proportional to the temporal 
and spatial distances between two points along the p3 and pi eikonals, and are conjugate to the energy and spatial 
momentum of emitted radiation. In the region that dominates the integrals, the conjugate of the spatial momenta 
is constrained to be smaller than the conjugate of the energy by a factor of /3, corresponding to a range in spatial 
momentum larger than energy by a factor of 1//3. Thus, when the energy of gluon exchange approaches the scale of 
heavy quark kinetic energy, m/3^, the exchanged momentum approaches m/3, which is the same order as the heavy 
quark momentum. In this region, the quarks can no longer by considered as recoilless sources, and the eikonal 
approximation is not reliable. In Ref. [7||, an effective theory treatment for resummation in this region was developed, 
and in ^Sj an equivalent analysis in which the pair is replaced by a single Wilson line was employed. In either case, it 
is necessary to assume that radiation with energy-momentum scales up to r7i/3 cancels in the sum over final states for 
the inclusive cross section. 

Finally, it is interesting to note that the connection between cos 9 and the kinematic variables s and ti used in the 
derivation of Eq. ((25|) is singular at absolute threshold. In two-to-two scattering processes, s + ti + ui =0, and the 



function Tg ' in Eq. ([M)) is determined by two variables, which we can choose as (s,ti). The absolute threshold 
/3 = then represents the single point (s = sq = Am^,ti = to = —sq/2) in the physical region of the kinematical 
(s,ti) plane; see Fig. [5] The relationship between ti and s can now be written as 
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FIG. 6: Illustration of the physical region (shaded) in the (s,ti) kinematic plane. The normalization is chosen such that so — 1. 
The solid lines depict the lines P — (i.e. s = so) and t = to- The dashed line denotes the set of points ii = — s/2 where the 
function Fg , defined in Eq. (I24p . vanishes due to antisymmetry (note that this line also crosses the absolute threshold point 



to 



-so/2) 



u 



'-{I- 13 COS0) 



fo 
2 



Vs 



\/s~~Sq COS! 



--is- 
2 ^ 



So) 



(27) 



At threshold, all derivatives of ii with respect to s diverge unless cos 9 — 0. Therefore, in terms of figure [51 all curves 
that approach absolute threshold at fixed, non-zero cos 6* converge to the vertical tangent to the boundary of the 
physical region at the point ii = to, s ~ sq. At the same time, it is straightforward to verify that along any straight 
line defined by the equation 



ti =to + c{s — So) , wi = to - (1 + c){s — So) , c < oo 



(28) 



the function T 



(134) 



in Eq. 



vanishes at threshold and thus matches smoothly to the behavior of the function at 
ui = ti (8|. In a sense, as /3 — >■ 0, the cosine becomes ill-defined, because it no longer uniquely defines the Mandelstam 
variables. At the same time, for any fixed value of P, and in particular, for /3 — ^ 1, the full expression for the soft 
anomalous dimension as found in Ref. [6] applies without subtlety for gluon radiation in the range identified above. 



V. CONCLUSIONS 

In this note we have demonstrated the mutual consistency of a full position-space calculation of the massive soft 
anomalous dimension matrix following Ref. [4|, and the momentum space results described in Ref. [6|. Inparticular, the 
color-symmetric structure of the two-loop double-exchange diagrams before renormalization shown in [J| is consistent 
with the results for planar diagrams found in Ref. |6|]) once renormalization is taken into account. We have also 
confirmed the consistency of results for the massive anomalous dimension matrix with those for the massless case and 
with the next-to-next-to leading logarithmic threshold resummations described in [7|, |8| . Finally, we have seen how 
coordinate space analysis provides a clear interpretation of the non-uniform limit of the anomalous dimension matrix 
at absolute threshold. 
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Appendix A: Explicit forms for the integrals 



The functions I{pi,pj,e), defined by Eq. ([3]), are one-dimensional integrals. For space-like kinematics (for example 
Pi incoming, pj outgoing) the result in d = 4 — 2e dimensions and for arbitrary masses can be expressed as 



I{Pi,P],£) 



r2i^^M2'^'2 



e;i 



9 9 

(pi -PjY 



(Al) 
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The above function has a weU defined massless hmit in 4 — 2^ dimensions: 2F1 (1/2, 1; 3/2 — e; 1) = —(1 — 2s)/{2e). 
As can be seen from Eq. (jAip all non-trivial dependence on the kinematical invariants and masses comes entirely 
through the combination: 



familiar from the one-loop case [2|. Expressions for the functions /(^^'^^ appearing in Eqs. pTj) . (IT51) can be obtained 
by expanding the hypergeometric function in e, or by expansion of the integrand in Eq. ([S]) before integration. For 
expansion of the full expression (|A1I) we have used the program HypExp [22| , as well as standard relations between 
the polylogarithmic functions. The resulting expressions can be written as 

The consistency of these results with those of Ref. 16] is readily checked by relating the variable /3.y = — cosh^^(±pi • 
Pjlniinij) to our Vij by 

A,.-iln(i^), coth/3.,.-l. (A5) 



While compact, the expressions in Eqs. (IA3|) and (jA4l) are somewhat inconvenient for analytic continuation from 
space-like to time-like kinematics, because their momentum-dependence is only through the squares [pi ■ pj)'^ in Vij. 
For completeness and for use in our discussion of the threshold limit, we recall that analytic continuation is made 
straightforward by re-expressing the amplitudes in terms of the variable 



(TTii+TTij)^ /-| {mi—mj)^ 






1 + Vij A {mj+mjy ^ L {m.-m,)^ 1 -|- xfj 



(A6) 



with Sij — {pi + Pj) ■ The case of rrii — nij — m is of particular interest, where 



1 4m^ 1 

i — i 1 L 

Sij 1 U'lj 



1- 



iml^l l + bij ' s^j (l-a;y)2' 



(A7) 



where bij is the center-of-mass velocity for the pair pi, pj. For arbitrary masses, we have Xij > for Sy < (space-like) 
and Xij < for Sij > (time- like). 

For any masses, we now rewrite the functions I^^^-^> in terms of the Xij as 

/(-i)(p„Pj) = --^^ ln(x,,), (A8) 

I^"HP^,P,) = \^ (-Li2(4) + ln2(x.,) - 21n(a;,,) ln(l - 4) + ^) . (A9) 

For space-like kinematics, where 1 > Xij > 0, these functions are real. Analytic continuation to time-like kinematics 
is found by the replacement: Xij — > — 1x^1 + ie = \xij\e^'" (see Sec. 6 of Ref. [2^ for details). Following this rule, the 
following terms should be added to Eqs. (jASp . (|A4p when the kinematics is time-like: 



A/(-^)(p„ft) = -^n\±^ = -^, (AlO) 

-L X-- Vj/i 



„2 



A/(o)(p„p,) - ^ii^ (-^2^i7r [21n|a;,,|-21n(l-4)] 



TT^ in 



Vi 



ln(i^^U21n^^ + "^^ 



1 -I- Vij ) \ 2vij 



(All) 
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Finally, for space-like kinematics the result for the non-planar diagram in Fig. [2l as expressed in Eq. (jTS]) above, 
reads [y]: 

/3s(k,P,) = -\n\x^,,) = -iln^ (^-p^) . (A12) 

In time-like kinematics, it also receives the following term: 

^ha{p^.P,) - 7r2 - 2i7rln(|xy I) = tt^ - ivrln {\-^\ ■ (A13) 



[1] N. Kidonakis and G. Sterman, Nucl. Phys. B 505, 321 (1997) [arXiv:hep-ph/9705234]; N. Kidonakis, G. Oderda and 
G. Sterman, Nucl. Phys. B 531, 365 (1998) [arXiv:hep-ph/9803241]. 

S. Catani, S. Dittmaier and Z. Trocsanyi, Phys. Lett. B 500, 149 (2001) [arXiv:hep-ph/0011222]. 
N. Kidonakis, Phys. Rev. Lett. 102, 232003 (2009) [arXiv:0903.2561 [hep-ph]]. 
A. Mitov, G. Sterman and I. Sung, Phys. Rev. D 79, 094015 (2009) [arXiv:0903.3241 [hep-ph]]. 
T. Becher and M. Neubert, Phys. Rev. D 79, 125004 (2009) [arXiv:0904.1021 [hep-ph]]. 

A. Ferrogha, M. Neubert, B. D. Pecjak and L. L. Yang, arXiv:0907.4791 [hep-ph]; arXiv:0908.3676 [hep-ph]. 
M. Beneke, P. Falgari and C. Schwinn, arXiv:0907.1443 [hep-ph]. 

M. Czakon, A. Mitov and G. Sterman, Phys. Rev. D 80 (2009) 074017 [arXiv:0907.1790 [hep-ph]]. 

V. Ahrens, A. Ferrogha, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Lett. B 687, 331 (2010) [arXiv:0912.3375 [hep- 
ph]]; V. Ahrens, A. Ferrogha, M. Neubert, B. D. Pecjak and L. L. Yang, arXiv: 1003.5827 [hep-ph]; N. Kidonakis, Phys. 
Rev. D 81, 054028 (2010) [arXiv:1001.5034 [hep-ph]]. 

L. J. Dixon, L. Magnea and G. Sterman, JHEP 0808, 022 (2008) [arXiv:0805.3515 [hep-ph]]. 

T. Becher and M. Neubert, Phys. Rev. Lett. 102 (2009) 162001 [arXiv:0901.0722 [hep-ph]]; JHEP 0906, 081 (2009) 
[arXiv:0903.1126 [hep-ph]]. 

E. Gardi and L. Magnea, JHEP 0903, 079 (2009) [arXiv:0901.1091 [hep-ph]]. 

S. M. Aybat, L. J. Dixon and G. Sterman, Phys. Rev. D 74, 074004 (2006) [arXiv:hep-ph/0607309] . 
G. Sterman and M. E. Tejeda-Yeomans, Phys. Lett. B 552, 48 (2003) [arXiv:hep-ph/0210130]. 

A. Mitov and S. Moch, JHEP 0705, 001 (2007) [arXiv:hep-ph/0612149]; T. Becher and K. Melnikov, JHEP 0706, 084 
(2007) [arXiv:0704.3582 [hep-ph]]. 

J. Gluza, A. Mitov, S. Moch and T. Riemann, JHEP 0907, 001 (2009) [arXiv:0905.1137 [hep-ph]]. 

S. Catani and M. H. Seymour, Phys. Lett. B 378, 287 (1996) [hep-ph/9602277]; Nucl. Phys. B 485, 291 (1997) [Err .-ibid. 
B 510, 503 (1997)] [hep-ph/9605323]. 

I. A. Korchemskaya and G. P. Korchemsky, Phys. Lett. B 287, 169 (1992); G. P. Korchemsky and G. Marchesini, Nucl. 
Phys. B 406, 225 (1993) [arXiv:hep-ph/9210281]. 
S. Catani, Phys. Lett. B 427, 161 (1998) [arXiv:hep-ph/9802439]. 
A. Sen, Phys. Rev. D 28, 860 (1983). 

K. Hagiwara, Y. Sumino and H. Yokoya, Phys. Lett. B 666, 71 (2008) [arXiv:0804.1014 [hep-ph]]; Y. Kiyo, J. H. Kuhn, 
S. Moch, M. Steinhauser and P. Uwer, Eur. Phys. J. C 60, 375 (2009) [arXiv;0812.0919 [hep-ph]]; M. Beneke, M. Czakon, 
P. Falgari, A. Mitov and C. Schwinn, arXiv:0911.5166 [hep-ph]. 

T. Huber and D. Maitre, Comput. Phys. Commun. 175, 122 (2006) [arXiv:hep-ph/0507094]; 
T. Huber and D. Maitre, Comput. Phys. Commun. 178, 755 (2008) [arXiv:0708.2443 [hep-ph]]. 

W. Bernreuther, R. Bonciani, T. Gehrmann, R. Heinesch, T. Leineweber, P. Mastrolia and E. Remiddi, Nucl. Phys. B 
706, 245 (2005) [arXiv:hep-ph/0406046]. 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 



[10 
[11 

[12; 

[13; 

[14 
[15 

[16 

[17 

[18 

[19 
[20 
[21 



[22; 
[23; 



